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Enhancements  to  the  early  version  of  the  COVER  and  BROWSE  programs  are  described  herein. 
In  addition,  there  are  changes  to  several  of  the  keyword  names  and  formats  linking  the  COVER 
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COVER 


SHRUB 


Invoke  the  COVER  option  in 
the  shrub  and  cover  extension, 
field  1:  Method  to  be  used 
to  compute  foliage  biomass. 
Invoke  the  BROWSE  option  of 
the  shrub  and  cover  extension, 
field  1:  Number  of  years 
since  stand  was  regenerated, 
field  2:  Number  of  years 
shrub  output  will  be  printed, 
field  3:  Habitat  type  code  for 
processing  SHRUB  option. 
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Keyword 

COVER         Invoke  the  COVER  extension 
(for  either  canopy  or  shrubs), 
field  1:  Cycle  to  begin 
COVER  predictions, 
field  2:  Dataset  reference 
number  for  output. 

CANOPY       Compute  crown  cover  statistics. 


SHRUBS       Compute  shrub  statistics. 

field  1:  Number  of  years 
since  stand  disturbance, 
field  2:  Habitat  type  code  for 
processing  SHRUBS  option, 
field  3:  Physiographic  type 
code. 

field  4:  Disturbance  type  code. 
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RESEARCH  SUMMARY 

The  COVER  extension  to  the  Stand  Prognosis  Model 
predicts  tree  canopy  closure,  crown  volunne,  crown  profile 
area,  and  foliage  biomass  within  vertical  height  classes, 
and  the  probability  of  occurrence,  height,  and  cover  of 
shrubs  in  forest  stands.  The  model  may  be  used  to  pro- 
duce a  descriptive  summary  of  a  stand  at  the  time  of 
inventory,  or  to  project  overstory  and  understory  charac- 
teristics through  time  for  natural  and  managed  stands. 
This  paper  documents  use  of  the  COVER  program,  an 
adjunct  to  Version  5.0  of  the  Prognosis  Model.  Prepara- 
tion of  input,  interpretation  of  output,  program  control, 
and  model  characteristics  are  described.  Potential  appli- 
cations of  COVER  estimates  to  wildlife,  hydrology,  and 
insect  pest  modeling  are  presented. 
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the  CANOPY  and  SHRUBS 
Extension  of  the  Stand 
Prognosis  IVIodel 

Melinda  Moeur 


INTRODUCTION 

Forest  managers  in  the  Northern  Rocky  Mountains  have  extensively  used  the 
Prognosis  Model  for  Stand  Development  (Stage  1973)  to  summarize  current  stand 
conditions,  and  to  predict  the  future  path  of  stand  growth  and  the  likely  conse- 
quences of  alternative  management  practices  on  stand  development.  The  use  of  stand 
models  like  Prognosis  need  not  be  restricted  to  timber  management  applications.  For 
example,  forest  managers  must  also  consider  how  management  practices  may  change 
the  suitability  of  a  stand  for  wildlife  habitat,  the  composition  of  the  understory,  or 
the  sequence  of  succession.  The  computer  model  described  in  this  paper,  known  as 
COVER,  extends  Prognosis  by  modeling  the  development  of  tree  crowns  and  under- 
story vegetation.  COVER  provides  three  types  of  information:  a  description  of  the 
amount  of  cover  and  foliage  in  the  tree  canopy  by  height  class;  the  height  and  cover 
of  shrubs,  forbs,  and  grasses  in  the  understory;  and  a  summary  of  overstory  and 
understory  cover  and  biomass  for  the  stand. 

COVER  offers  two  options  that  can  be  run  separately  or  together:  CANOPY  con- 
trols predictions  of  values  related  to  tree  crowns;  SHRUBS  controls  predictions  of 
understory  characteristics.  Possible  applications  of  the  COVER  model  include 
examining  the  likely  effects  of  silvicviltural  treatments  on: 

•  forest  stand  characteristics  important  to  wildlife,  such  as  thermal  cover  and  hiding 
cover,  browse  production,  and  the  interactions  of  shrubs  and  trees  that  determine 
vertical  and  horizontal  stand  structure. 

•  dynamics  of  the  shrub  community  affecting  stand  succession  and  competition  with 
regeneration. 

•  vertical  crown  form  and  foliage  distribution  important  to  feeding  and  dispersal 
patterns  of  insect  pests,  such  as  the  western  spruce  budworm. 

•  canopy  cover  and  ground  cover  development  which  affect  the  hydrologic  charac- 
teristics of  a  stand  following  harvest. 

This  manual  is  intended  to  be  a  companion  publication  to  the  Prognosis  Model 
User's  Guide  (Wykoff  and  others  1982).  The  research  described  here  is  based  on  the 
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idea  that  management  prescriptions  which  affect  nonvegetation  ecosystem  compo- 
nents can  be  compared  and  evaluated  by  examining  simulated  changes  in  the  major 
vegetation  components.  Wykoff  and  others  (1982)  have  summed  up  this  philosophy, 
which  guides  our  research: 

Consequences  for  streamflow  from  the  forest,  for  wildlife  populations, 
and  for  pest  populations  that  inhabit  the  forest,  as  well  as  the  capabil- 
ity of  the  forest  to  yield  timber  or  provide  recreation— all  depend  on 
how  the  dominant  vegetation  changes  and  is  changed.  Unfortunately, 
yield  forecasts  have  traditionally  emphasized  the  merchantable  harvest 
that  might  be  obtained,  either  immediately  or  as  a  sequence  of  yields 
obtainable  at  intervals  of  time  into  the  future.  Volumes  of  merchant- 
able timber  have  been  the  most  common  units  of  measure  because 
timber  products  have  usually  been  the  primary  reason  for  investment. 
As  other  uses  for  the  forest  become  more  important,  however,  growth 
forecasts  need  to  be  stated  in  more  fundamental  descriptions  of  the 
future  forest  stand.  Too  often,  evaluation  of  tradeoffs  among  conflict- 
ing activities  or  objectives  for  use  of  forest  resources  has  been  ham- 
pered by  lack  of  sensitivity  of  the  forecasts  to  the  interactions  among 
ecosystem  components.  One  objective  for  development  of  the  Stand 
Prognosis  Model  is  to  so  characterize  stand  dynamics  that  the  model 
will  provide  a  sensitive  basis  for  representing  interactions  involving 
the  tree  species. 

COVER  can  provide  a  detailed  picture  of  the  vegetative  structure  of  a  single  stand 
through  time  under  different  sUvicultural  prescriptions.  By  linking  to  the  Parallel 
Processing  version  of  Prognosis  (Crookston  1985),  COVER  can  also  model  long-term, 
large-scale  changes  for  groups  of  stands  arranged  in  time  and  space.  As  such,  it  may 
be  a  useful  tool  for  enhancing  forest  management  decisions  that  concern  nontimber 
ecosystem  components.  This  publication  begins  with  a  description  of  the  information 
produced  by  the  model  and  instructions  for  making  the  program  run.  There  follows  a 
discussion  of  the  biological  behavior  of  the  individual  submodels  and  a  final  section 
dealing  with  potential  applications. 

Range  of  Predictions         The  models  comprising  COVER  are  parameterized  with  data  collected  in  the 

Inland  Northwest  and  Northern  Rocky  Mountain  forests.  The  user  should  determine 
if  the  range  of  species  and  site  conditions  for  which  predictions  are  made  are 
applicable  to  the  local  situation. 

Submodels  for  the  CANOPY  option  predict  conifer  crown  width,  crown  shape, 
and  fohage  biomass.  Equations  for  crown  width  and  foliage  biomass  are  derived 
from  data  on  370  trees  on  14  sites  in  northern  Idaho  and  western  Montana 
(Moeur  1981;  Brown  1978).  Sampled  stand  basal  area  ranges  from  1  to 
426  ft^/acre.  Open-grown  trees  and  trees  that  were  obviously  damaged  or 
heavily  defohated  were  not  sampled.  Data  for  the  crown  shape  models  are  from 
9,800  trees  on  12  sites  in  western  Montana,  the  Blue  Mountains  of  eastern 
Oregon,  and  the  University  of  Idaho  forest  near  Moscow,  ID  (Langelier  and 
Garton  in  press  a). 

The  understory  portion  of  COVER  includes  models  that  predict  probability  of 
occurrence,  height,  and  cover  of  individual  shrub,  forb,  grass,  and  fern  species 
(Scharosch  1984;  Laursen  1984).  Understory  data  are  from  over  10,000 
1/300-acre  plots  on  500  stands  in  Douglas-fir,  grand  fir,  western  redcedar,  west- 
ern hemlock,  and  subalpine  fir  habitat  types.  Stands  were  measured  between  3 
and  40  years  following  major  stand  disturbance.  These  data  are  from  northern 
Idaho,  northeastern  Washington,  and  northwestern  Montana  (Colville, 


2 


Panhandle,  Kootenai,  Lolo,  and  Clearwater  National  Forests),  to  central  and 
southern  Idaho  and  northwestern  Wyoming  (Nezperce,  Boise,  Payette,  and 
Targhee  National  Forests)  (Ferguson  and  others  [in  press]).  In  addition,  Irwin 
and  Peek  (1979)  fit  models  for  twig  production  and  dormant  season  shrub 
biomass  on  a  subset  of  the  data  (grand  fir,  cedar,  and  hemlock  types). 

Current  data  sources  make  the  COVER  predictions  most  applicable  to  the 
Northern  Region  (R-1),  and  portions  of  the  Intermountain  (R-4)  and  Pacific 
Northwest  (R-6)  Regions  of  the  Forest  Service.  The  user  should  exercise  caution 
in  extending  predictions  outside  these  geographic  areas. 

Data  Requirements  Information  needed  to  run  COVER  consists  of  the  minimum  Prognosis  Model 

input— the  inventory  design  used  to  measure  the  stand,  a  list  of  sampled  trees  for 
which  species,  diameter,  and  plot  identification  have  been  recorded,  and  values  for 
slope,  aspect,  elevation,  habitat  type,  and  forest  location  recorded  on  the  STDINFO 
card  (Wykoff  and  others  1982).  If  the  SHRUBS  option  is  used,  time  since  stand 
disturbance,  type  of  disturbance,  and  physiographic  position  are  reqiiired. 

Understory  predictions  are  improved  if  field  measurements  of  shrub  height  and 
cover  are  available  for  calibrating  portions  of  the  shrub  models.  These  data  are  only 
supplementary,  and  both  the  CANOPY  and  SHRUBS  options  will  execute  without 
them.  The  keywords  section  discusses  how  to  enter  calibration  information. 

INFORMATION  PRODUCED 

The  COVER  program  normally  produces  three  displays— one  describing  the 
structure  of  tree  crowns,  another  describing  the  composition  of  the  understory, 
and  a  third  summarizing  overstory  and  understory  cover  and  biomass.  The  user 
may  insert  keywords  to  turn  off  the  printing  of  any  of  the  displays. 

As  you  proceed  through  this  publication,  the  same  stand  used  in  the 
Prognosis  Model  User's  Guide  (Wykoff  and  others  1982)  vdll  be  used  to  develop 
examples.  The  stand  (S248112)  is  on  the  St.  Joe  National  Forest.  It  is  57  years 
old  at  the  inventory  date,  positioned  on  a  northwest  aspect,  30  percent 
midslope,  at  3,400  feet  of  elevation,  and  is  a  Tsuga  heterophylla/Clintonia 
uniflora  habitat  type.  Values  presented  here  result  from  the  COVER  extension 
combined  with  Version  5.0  of  the  Inland  Empire  Prognosis  Model  and  Version 
1.0  of  the  Regeneration  Estabhshment  Model  (Ferguson  and  Crookston  1984). 
Four  COVER  extension  keywords  (COVER,  CANOPY,  SHRUBS,  and  END) 
inserted  into  the  Prognosis  Model  runstream  in  figure  1  produce  the  example 
canopy  and  shrub  displays. 

The  Canopy  Cover  Statistics  display  is  the  first  of  the  three  COVER  displays 
(fig.  2).  It  is  produced  when  the  CANOPY  option  is  specified.  To  give  the  user  a 
feel  for  the  vertical  profile  of  the  conifer  component  of  the  stand,  crown  cover 
values  are  partitioned  by  10-foot  height  classes. 

Trees  per  acre.— The  number  of  trees  per  acre  whose  total  heights  fall  within 
a  given  height  class. 

Canopy  closure.— The  percentage  of  ground  area  covered  by  the  projections  of 
individual  crowns  of  trees  whose  total  heights  fall  within  a  given  height  class 
(fig.  3a). 

Crown  profile  area.— The  area  in  square  feet  per  acre  within  vertical  height 
classes  occupied  by  crown  profiles,  represented  by  the  sum  of  lateral  areas  of 
crown  profile  sections  within  height  classes  (fig.  3b).  Crown  profile  area  may  be 
thought  of  as  a  foliage-height  profile,  or  the  view  one  would  have  if  the  stand 
were  "squashed  flat"  in  a  vertical  plane. 


Canopy  Cover 
Statistics  Display 


3 


USER'S  MANUAL  EXAMPLE  STAND 
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Figure  1  .—Keyword  and  tree  record  file  used  to  project  stand  S248112 
using  the  combined  Prognosis  Modei  and  COVER  extension. 
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Crown  volume.— The  volume  in  cubic  feet  per  acre  of  tree  crowns  within 
height  classes.  Volume  is  determined  using  integration  of  standard  volume 
formulas  for  conical,  spherical,  parabolic,  elliptic,  and  neiloid  crown  forms. 

Foliage  biomass.— The  biomass  of  foliage  in  pounds  per  acre  in  the  stand, 
represented  by  the  sum  of  foliage  in  individual  crown  sections  within  height 
classes. 

Shrub  Statistics  The  second  COVER  display,  shrub  statistics  by  species  (fig.  4)  is  produced 

Display  when  the  SHRUBS  keyword  is  present  in  the  runstream.  Predictions  are  made 

for  each  of  31  species  or  species  groups  hsted  in  table  1.  Nine  species  with  the 
greatest  predicted  cover  in  the  stand  are  always  displayed:  three  each  within 
low  (less  than  1.7  ft),  medium  (1.7  —  7  ft),  and  tall  (7  ft  and  greater)  height 
classes  (Patterson  and  others  1985). 

Cover.— The  percentage  of  area  covered  by  the  canopy  of  an  individual 
species  projected  vertically  onto  the  ground.  Species  cover  is  weighted  by 
probability  of  occurrence. 
Height.— Average  height  of  the  species,  in  feet. 

Probability  of  occurrence.— The  proportion  (expressed  as  a  percentage)  of 
1/300-acre  plots  in  the  described  stand  that  contain  the  indicated  species. 

The  remaining  species  within  each  height  class  are  combined  into  an  "other" 
category.  For  these,  cover  is  the  sum  of  other  species  cover  weighted  by  their 
probabilities,  height  is  their  average  height  weighted  by  species  cover, 
and  probability  is  the  sum  of  their  individual  probabilities.  In  addition,  the  user 
may  specify  up  to  six  species  which  will  always  be  displayed,  using  the 
SHOWSHRB  keyword  card  and  supplemental  data  record. 

Shrub  values  are  computed  only  between  3  and  40  years  following  stand 
disturbance  (the  lower  and  upper  limits  of  the  range  of  data  used  to  model 
shrub  production).  In  the  course  of  a  projection,  if  a  thinning  causes  canopy 
closure  to  fall  below  50  percent,  and  if  the  thinned  volume  is  at  least  20  percent 
of  the  volume  before  thinning,  then  time  since  stand  disturbance  is  reset  to  3 
years  and  shrub  computations  resume.  In  the  example  simulation,  time  since 
disturbance  was  incremented  from  initial  stand  age  of  57  years  until  the 
overstory  removal  in  year  2037.  Neither  the  thinning  in  1977  nor  the  one  in 
2007  opened  the  canopy  enough  to  trigger  shrub  calculations  in  the  program. 
The  overwood  cut  in  2037  reduced  crown  closure  from  98  percent  to  14  percent, 
causing  the  shrub  response  calculations  to  begin. 
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STAND  GROWTH  PROGNOSIS  SYSTEM  VERSION  5.0  --    INLAND  EMPIRE 

STAND   ID:     52146112  MANAGEMENT  CODE:     NONE  PROGNOSIS  WITH  COVER  EXTENSION  -  USER'S  MANUAL  EXAMPLE  STAND 


SHRUB  STATISTICS 


LOW  SPEC  I ES   (0-1.7   FT  ) 


ARUV: ARCTOS lAPHYLOS  UVA-URS I 
BERB: BERBER  I S  SPP. 
LIBO:LINNAEA  BOREAL  I S 
PAMY: PACH I  ST  I  MA  MYRSINITES 
SPBE:SPIRAEA  BETULIFOLIA 
VASC: VACCI N I UM  SCOPARIUM 
CARX:CAREX  SPP. 


MED  I UM  SPEC  I ES   (1.7-7   FT  ) 


LON I : LON I  CERA  SPP. 
MEFE: MENZ I ES I  A  FERRUGINEA 
PHMA: PHYSOCARPUS  MALVACEUS 
R I  BE: R I BES  SPP. 
ROSA:ROSA  SPP. 
RUPA:RUBUS  PARVIFLORUS 
SMCA: SHEPHERD  I  A  CANADENSIS 
SYMP: SYMPHOR I CARPOS  SPP. 
VAME: VACCI N I UM  MEMBRANACEUM 
XETE:XEROPHYLLUM  TENAX 
FERN: FERNS 

COMB:OTHER  SHRUBS  COMBINED 


TALL   SPECIES   (7+  FT) 


ACGL:ACER  GLABRUM 
ALSI :ALNUS  SINUATA 
AMAL:AMELANCHI ER  ALNIFOLIA 
CESA:CEANOTHUS  SANGUINEUS 
CEVE:CEANOTHUS  VELUTINUS 
COST:CORNUS  STOLON  I FERA 
HOD  I : HOLOD I SCUS  DISCOLOR 
PREM:PRUNUS   EMARG I NATA 
PRV I : PRUNUS  V  I RG I N I  ANA 
SALX:SALIX  SPP. 
SAMB: SAMBUCUS  SPP. 
SORB:SORBUS  SPP. 


ATTRIBUTES  OF  THE  FIRST   THREE  SPECIES  WITH  GREATEST  COVER   IN  EACH  HEIGHT  GROUP 
(ALL  OTHERS  WITHIN  GROUP  COMBINED   INTO  CATEGORY  "OTHR"  ) 


COVER  --  SPECIES  COVER 
HEIGHT  --  AVERAGE  SPECIES  HEIGHT  (FEET) 

PROB  --  SPECIES  PROBABILITY  OF  OCCURRENCE 


YEAR 

LOW  SPECIES 

MED  1 UM 

SPECI ES 

TALL 

SPEC 

ES 

1977 

T  1  ME 

S  1  NCE 

Dl STURBANCE= 

57 

EXCEEDS 

14  0 

YEARS 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

1977  POST- 

TH  1  N 

T  1  ME 

S  1  NCE 

Dl STURBANCE= 

57 

EXCEEDS 

'(0 

YEARS. 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

1987 

TIME 

S  1  NCE 

D  1  STURBANCE-= 

67 

EXCEEDS 

140 

YEARS. 

SHRUB 

STATISTICS 

NOT 

COMPUTED 

1997 

T  1  ME 

SI  NCE 

D 1 STURBANCE= 

77 

EXCEEDS 

40 

YEARS. 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

2007 

TIME 

SINCE 

D 1 STURBANCE= 

87 

EXCEEDS 

140 

YEARS. 

SHRUB 

STATISTICS 

NOT 

COMPUTED 

2007  POST- 

TH  1  N 

T  1  ME 

SINGE 

D 1 STURBANCE= 

87 

EXCEEDS 

140 

YEARS. 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

2017 

T  IME 

SINCE 

Dl STURBANCE= 

97 

EXCEEDS 

140 

YEARS. 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

2027 

T  IME 

S  1  NCE 

D 1 STURBANCE= 

107 

EXCEEDS 

140 

YEARS. 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

2037 

T  1  ME 

S  1  NCE 

D 1 STURBANCE= 

117 

EXCEEDS 

140 

YEARS. 

SHRUB 

STAT  1  ST  1 CS 

NOT 

COMPUTED 

2037: POST- 

TH  1  N 

SPEC  1  ES 

LI  BO     SPBE      PAMY  OTHR 

VAME  FERN 

RUPA 

OTHR 

ACGL  COST 

AMAL  OTHR 

COVER 

7 

5  2 

3       0.7  0. 

0 

13.2  t| 

9 

14.  3 

10.5 

11.7  0.8 

0 

6  0.8 

HEIGHT 

0. 

5  1 

8       1.5  0. 

7 

2  .  U  2 

1 

2.3 

2.9 

8.3  5.5 

5 

5  4.7 

PROB 

27. 

7  9. 

6       5.1  0. 

0 

50.3  25 

2 

26.9 

71.7 

39.0  2.2 

14 

2  8.9 

2047 
SPEC  I ES 
COVER 
HE  I GHT 
PROB 


L  I  BO 
8.9 
0.5 

31.3 


PAMY 
5.8 
1  .  7 

314.  2 


SPBE 
2.2 
1  .8 
8.6 


OTHR 
0.  1 
0.7 
O.I4 


VAME  RUPA 

8.7  I4.  1 

2.I4  2.14 

37.9  22.6 


LON  I 
2.5 
3.  1 

16.14 


OTHR 
7.5 
2.8 

I42.  8 


ACGL 
14. 9 
8.8 

15.5 


SALX 
1  .0 
8.9 
5.2 


AMAL 
0.9 
6.6 
5.9 


OTHR 
2.0 
5.9 

11.2 


2057 
SPEC  I ES 
COVER 
HEIGHT 
PROB 


LI  BO 
9.9 
0.5 

32.9 


PAMY 
7.3 
1  .6 

141  .  14 


SPBE 
2.  1 
1  .8 
7.9 


OTHR 
0.  1 
0.7 
0.6 


VAME 
6.6 
2.3 

33.9 


RUPA 
3  .  1 
2.3 

18.3 


LON  I 
2.  1 
3.0 

114.  5 


OTHR 
5.7 
2.8 

3I4.6 


ACGL 
3.0 
8.3 

10.9 


SALX 
1  .  3 
9.2 
7.7 


AMAL 
0.9 
6.5 
5.8 


OTHR 
1  .5 
6.1 

9.3 


2067 
SPECI ES 
COVER 
HE IGHT 
PROB 

20  7  7 
SPEC  I ES 
COVER 
HE IGHT 
PROB 


L  I  BO 
10.7 
0.5 
33  . 14 


L  I  BO 
11.5 
0.5 
33 


PAMY 
6.8 
1  .  5 

39.2 


PAMY 
5.7 
1  . 14 

3I4.2 


SPBE 
1  .  8 
1  .  7 
6.9 


SPBE 
1  .5 
1  .6 
5.9 


OTHR 
0  .  2 
0.7 
0.8 


OTHR 
0.2 
0.  7 
1  .  1 


VAME  RUPA  LON I  OTHR 

14.8  2.1  1.5  14.  3 

2.3  2.2  2.8  2.7 

30.3  13.7  11.8  28.8 


VAME 
3. 14 
2.2 
26.8 


RUPA 
1  .  2 
2.0 
9.6 


ROSA 
1  .  1 
2.8 
10.1 


OTHR 
3.  1 
2.5 
23.0 


ACGL 
1  .8 
7.3 
8.  1 


ACGL 
1  .  0 
6.3 
6 


SALX  AMAL 

1.2  0.7 

8.7  5.7 

8.0  5.0 


SALX 
1  .0 
7.7 
7.6 


AMAL 
0.5 
14. 8 
14. 0 


^,     PROB  33.9     314.2      5.9       1.1         26.8      9.6     10.1     23.0  6.0       7.6  4.0 

Figure  4.  — The  shrub  statistics  display  is  produced  when  the  SI-IRUBS  l<eyword  is  present  in  the  runstream. 


OTHR 
0.9 
6.0 
6.7 


OTHR 
0.5 
5.9 
4.7 


9 


Table  1.— Understory  species  for  which  predictions  are  currently  made  in  the  SHRUBS  portion  of 
the  COVER  program,  height  class,  and  source  of  information 


Code 

Scientific  name 

Common  name 

Height  class^ 

Source^ 

ACGL 

Arpf  n l^hn im 

Rnrkv  Mniintain  manip 

J 

a  b 

ALSI 

1 1  I  {J  J    \j  1 1  1  U  Ct  I  CI 

T 

1 

a 

AMAL 

oci  viucuci  1  y 

T 
1 

a,  U 

ARUV 

Arctostaphylos  uva-ursi 

Kinnikinnick 

L 

a 

SERB 

\_/  1  ^  ^      1  1    ^  1  ci  ^ 

l_ 

CI 

CARX 

r^a  rp  y  q  n  n 

L 

d 

CESA 

<w/  c  CJ 1  (t^if  1  KJ  J   •ja  1  1  ^  u  f  1  )  cr  u  o 

Rpdc;tpm  pPAnnthii'^ 

J 

n  h 

CEVE 

Shinvlpaf  rPAnnthii^^ 

Ol  Ml  ly  I^Cll    ^COI  IwLI  lUO 

J 

CI,  u 

COST 

J 

a 

FERN 

A  th\/ri  1  im  filiy-fpminj^ 
Ptsridi ufTi  BQuilinum 

FGrn 

M 

HODI 

Holodiscus  discolor 

Ocesn-spray 

J 

a,  b 

LIBO 

T\A/i  n  f  1  nvA/p  r 

1  VV  II  1  1  1 W  vV  ^ 1 

|_ 

LONI 

1  nnirpr^  ^nn 

H  nnpvc;  1  IP  k  lp 

1  iwii^yoLJwrxio 

M 

a.  b 

MEFE 

ivi  a          1  a  Id  1  L/y/  /  icct 

M  p  n  7  i  p  Q  i  3 

IVI  CI  1^  1  Co  1  CI 

M 

a 
d 

RAM  Y 

iCI^IIIOlllllCl 

[_ 

PH  MA 

iliyjKj\jaffJ\J^   1 1 1  a  1  V  Ci^  C/  U  J 

M  i  nph?i  rk 

1  N  1  1  1  C     Cl  1  r\ 

M 

A  h 

PREM 

1    1  {J  1 1  Li  J  C^IIICil^lllCilCl 

RittprrhPTfu 
Lji  iici      Id  1  y 

J 

CI,  u 

PRVI 

P run  US  vitQiniana 

Common  chokechGrry 

J 

RIBE 

f?/'Hpe  cnn 

0\  1  rra  nt 
w  u 1 1  di  1 1 

M 

a  b 

ROSA 

Rn<^^  cnn 
n  L/O  C3   o\J\J . 

Ross 

M 

RUPA 

Rubus  parviflorus 

Thimbleberry 

M 

a,  b 

SALX 

Salix  spp. 

Willow 

T 

a!  b 

SAMB 

Sambucus  spp. 

EldsrbGrry 

T 

SHCA 

Rii^^ptt  hiiffalnhpcrv 

M 

SORB 

w  \^  1     \j  \j   <j  . 

Mou  ntai  n-3sh 

T 

a 

SPBE 

Spiraea  betulifolia 

Shinyleaf  spiraea 

L 

a.  b 

SYMP 

Symphoricarpos  spp. 

Snowberry 

M 

a,  b 

VAME 

Vaccinium  membranaceum 
Vaccinium  globulare 

Big  huckleberry 
Globe  huckleberry 

M 

a.  b 

VASC 

Vaccinium  scoparium 

Grouse  whortleberry 

L 

a 

XETE 

Xerophyllum  tenax 

Common  beargrass 

M 

a 

COMB 

Artemisia  tridentata 
Clematis  columbiana 
Cornus  nuttallii 
Crataegus  douglasii 
Juniperus  spp. 
Ledum  glanduiosum 
Lonicera  caerulea 
Lonicera  in  value  rat  a 
Opiopanax  horridum 
Ptiiladelphus  lewisii 

Other  shrubs  combined 
Prunus  pensylvanica 
Purshia  tridentata 
Rhamnus  purstiiana 
Rhododendron  albiflorum 
Rhus  trilobata 
Rubus  leucodermis 
Rubus  ursinus 
Spiraea  pyramidata 
Tax  us  brevi  folia 
Vaccinium  caespitosum 

M 

a 

''T  =  tall,  M  =  medium,  L  =  low  (from  Patterson  and  others  1985). 
^a:  Probability  of  occurrence  from  Scharoscfi  (1984). 

Height  and  percent  cover  from  Laursen  (1984). 
b:  Twig  production  and  dormant  season  biomass  (ABGR/CLUN,  THPL/CLUN,  TSHE/CLUN 

habitat  types  only)  from  Irwin  and  Peel^  (1979). 
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Canopy  and  Shrubs         The  third  display  is  the  Canopy  and  Shrubs  Summary  display  (fig.  5). 

Summary  Display  Understory  attributes: 

Time  since  disturbance.— Time  in  years  since  the  stand  has  been  entered  for 
harvest  or  site  preparation.  At  the  start  of  a  projection,  this  value  is  set  equal 
either  to  stand  age  or  to  the  value  entered  on  the  SHRUBS  keyword  record, 
and  is  incremented  by  the  length  of  each  cycle.  A  thinning  which  reduces 
canopy  closure  below  50  percent  and  accounts  for  at  least  20  percent  of  the 
prethinning  volume  causes  time  since  disturbance  to  be  reset  to  3  years.  See 
the  previous  section  for  more  explanation. 

Probability  of  shrub  cover  being  greater  than  zero.— The  proportion  (ex- 
pressed as  a  percentage)  of  1/300-acre  plots  in  the  described  stand  that  contain 
any  shrub  cover. 

Shrub  cover.— Sum  of  individual  percentage  cover,  weighted  by  probability, 
for  all  species  in  low  (0  —  1.7  ft),  medium  (1.7  —  7  ft),  and  tall  (over  7  ft)  height 
classes,  and  total  shrub  cover. 

Average  shrub  height.— Average  height  in  feet  of  all  species,  weighted  by 
predicted  species  cover. 

Dormant  shrub  biomass.— Total  dormant  season  aboveground  shrub  biomass 
in  pounds  per  acre. 

Twig  production.— Current  year's  number  of  twigs  per  square  foot. 
Shrub  biomass  and  twig  production  are  only  computed  for  ABGR/CLUN, 
THPL/CLUN,  and  TSHE/CLUN  habitat  types  (codes  520,  530,  and  570). 

Stand  successional  stage  code.— A  vegetation  hfe  form  classification  described 
by  Peterson  (1982),  which  provides  a  basis  for  relating  wildlife  use  to  the  shrub 
and  conifer  structure  of  the  stand. 

Overstory  attributes: 

Stand  age.— Overstory  age  in  years,  entered  on  the  STDINFO  card  and  in- 
cremented by  the  length  of  each  projection  cycle.  The  value  of  stand  age  will 
change  during  the  projection  if  RE  SET  AGE,  a  regeneration  estabhshment 
model  keyword  (Ferguson  and  Crookston  1984),  is  used. 

Top  height.— Current  average  height  in  feet  of  the  largest  40  trees  per  acre 
by  d.b.h. 

Canopy  closure.— Total  percentage  canopy  closure. 

Crown  foliage  biomass.— Total  foliage  biomass  in  pounds  per  acre. 

Sum  of  stem  diameters.- The  sum  of  stem  diameters  at  breast  height  for  all 
trees  in  the  stand,  in  feet.  This  quantity  may  be  used  to  compute  stem  area 
available  for  hiding  cover  by  multiplying  by  height  value  of  interest  (for  exam- 
ple, average  shoulder  height  of  an  elk). 

Number  of  stems.— Total  number  of  trees  per  acre. 

Shrub-small  conifer  competition: 

If  the  SHRUBS  option  is  in  effect,  a  display  similar  in  format  to  the  Canopy 
Cover  table  is  produced.  It  expands  the  resolution  of  the  first  20  ft  of  the 
stand,  and  displays  both  shrub  cover  and  number  of  trees  cumulatively  by 
height. 

Shrub  cover.— Total  cover  of  shrubs  whose  predicted  heights  are  greater  than 
the  current  height. 

Number  of  trees.— Total  number  of  trees  per  acre  whose  heights  are  greater 
than  the  current  height. 
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PROGRAM  ORGANIZATION 

Each  projection  cycle  in  the  Prognosis  Model  starts  with  simulated  thinnings 
if  thinning  activities  have  been  scheduled  for  that  cycle.  If  either  the  CANOPY 
or  SHRUBS  option  has  been  selected,  crown  and  shrub  statistics  are  computed 
and  displayed  following  the  thinning.  Next,  diameter  and  height  growth,  change 
in  crown  ratio,  and  mortality  rate  are  computed  for  each  tree  record  in  the  inven- 
tory by  the  main  program.  The  tree  attributes  are  updated  at  the  end  of  the 
cycle  and  then  crown  and  shrub  statistics  are  computed  once  again.  Thus,  in 
thinning  cycles,  COVER  values  are  computed  and  displayed  twice.  In  projec- 
tion cycles  with  no  thinning,  there  is  only  one  call  to  the  COVER  extension  at 
the  end  of  the  cycle. 

Crown  cover  and  shrub  predictions  are  displayed  following  a  thinning  to  show 
the  immediate  effects  of  treatment  on  cover  development.  Predictions  and  dis- 
play for  crowns  and  for  shrubs  are  handled  somewhat  differently  in  the  post- 
thin  sequence.  In  the  CANOPY  subprogram,  the  crown  models  use  prethin 
levels  of  stand  density  statistics  (trees  per  acre,  basal  area,  and  relative  diam- 
eter) to  predict  crown  dimensions  on  remaining  trees.  This  is  done  because  no 
immediate  response  in  crown  dimension  or  foliage  on  residual  trees  to  a  change 
in  stand  density  would  be  expected.  In  the  SHRUBS  subprogram,  predictions 
are  made  using  postthin  stand  density  statistics  to  reflect  the  more  rapid  re- 
sponse of  understory  plants  that  one  would  expect  following  an  opening  up  of 
the  canopy. 

The  COVER  extension  may  be  used  in  conjunction  with  ESTAB,  a  Prognosis 
Model  extension  that  simulates  regeneration  and  subsequent  development  of  a 
regenerated  stand  (Ferguson  and  Crookston  1984).  Inside  the  combined 
COVER,  ESTAB,  and  Prognosis  Models,  shrub  values  are  computed  prior  to 
the  estabHshment  of  new  trees.  Although  presently  ESTAB  does  not  include  ef- 
fects of  shrub  development  on  the  establishment  of  new  trees,  the  computation 
sequence  has  been  set  up  in  anticipation  of  feedback  between  the  shrubs  and 
regeneration  systems.  For  now,  the  shrub-small  conifer  display  in  the  Summary 
table  may  help  the  user  examine  effects  of  shrub  competition  on  newly  estab- 
lished trees. 

Shrub  statistics  (probability  of  occurrence,  height,  cover,  biomass,  and  twig 
production)  are  not  computed  if  time  since  disturbance  exceeds  40  years.  If  the 
COVER  extension  is  used  in  conjunction  with  the  regeneration  establishment 
model,  and  stand  age  is  reset  to  a  value  less  than  40  years  using  the 
RESETAGE  keyword,  shrub  calculations  will  resume.  Also,  as  discussed  previ- 
ously, shrub  statistics  will  resume  beginning  with  the  postthin  predictions  for  a 
cycle  in  which  a  heavy  thinning  takes  place.  Canopy  cover  statistics  are  dis- 
played regardless  of  the  value  for  time  since  disturbance. 

KEYWORD  DESCRIPTIONS 

The  keyword  system  used  in  the  COVER  extension  is  similar  to  that  used  for 
the  Prognosis  Model.  Presently,  13  keywords  are  used  to  invoke  the  CANOPY 
or  SHRUBS  options,  supply  needed  information  to  the  extension,  or  modify  its 
output.  The  position  of  these  cards  in  the  Prognosis  Model  deck  is  unimportant 
except  that  they  should  be  in  a  group  beginning  with  the  COVER  keyword  and 
ending  with  the  END  keyword.  Like  aU  other  keywords,  they  must  precede  the 
PROCESS  card.  Formatting  of  keyword  records  follows  the  scheme  used  for  all 
other  Prognosis  Model  keywords.  Columns  1  to  10  are  reserved  for  the  keyword 
itself,  followed  by  seven  parameter  fields  of  10  columns  each.  Numeric  data 
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should  be  right-justified  within  the  parameter  field,  or  include  a  decimal  point. 
Three  of  the  keyword  records  in  the  SHRUBS  option  are  followed  by  supple- 
mental data  records.  The  END  record  signifies  the  end  of  keywords  for  the  ex- 
tension and  returns  control  to  the  main  program.  The  appendix  summarizes 
keyword  records  available  in  the  COVER  option. 

COVER  Keyword 

The  COVER  keyword  record  signifies  the  beginning  of  keywords  for  the 
extension.  It  has  two  parameter  fields. 

COVER     field  1:  Cycle  number  in  which  COVER  calculations 

begin;  default  =  beginning  of  projection.  COVER 
calculations  will  be  performed  in  all  cycles  subsequent  to  the 
specified  one. 

field  2:  Dataset  reference  number  for  COVER  output; 
default  =  18. 

CANOPY  Keyword 

The  CANOPY  keyword  invokes  the  crown  cover  options  of  the  extension, 
which  compute  crown  width,  shape,  and  foliage  biomass  for  each  tree  record, 
and  tree  cover  summary  statistics  for  the  stand.  It  has  no  parameter  fields 
associated  with  it. 

Understory  Options      SHRUBS  Keyword 

The  SHRUBS  keyword  tells  the  program  to  compute  shrub  statistics.  It  has 
four  parameter  fields. 

SHRUBS     field  1:  Time  in  years  since  stand  disturbance;  default  = 

stand  age.  If  stand  age  is  not  supphed,  default  =  3  years, 
field  2:  Habitat  type  code  selected  for  processing  shrub  options. 

Table  2  lists  habitat  types  for  which  shrub  predictions  are 
made;  default  =  stand  habitat  type  code, 
field  3:  Physiographic  type  code.  1  =  bottom,  2  =  lower  slope.  3  = 

midslope,  4  =  upper  slope,  5  =  ridge;  default  =  2. 
field  4:  Disturbance  type  code.  1  =  none,  2  =  mechanical,  3  = 
burn,  4  =  road;  default  =  1. 
As  indicated,  each  of  the  SHRUBS  keyword  parameters  has  a  default  value 
in  the  event  that  no  value  is  supplied  by  the  user.  The  value  entered  for  time 
since  disturbance  should  be  the  number  of  years  since  the  stand  was  entered 
for  harvest  or  site  preparation.  Time  since  disturbance  will  be  set  to  stand  age 
entered  on  the  STDINFO  keyword  record  if  the  user  fails  to  supply  a  value  on 
the  SHRUBS  keyword,  and  will  be  incremented  by  the  length  of  each  projec- 
tion cycle.  Although  stand  age  is  not  used  by  the  Prognosis  Model  to  calculate 
tree  growth,  time  since  disturbance  is  a  significant  predictor  of  shrub  develop- 
ment. As  noted  previously,  shrub  calculations  are  performed  only  if  the  current 
value  of  time  since  disturbance  is  between  3  and  40  years.  Three  situations 
may  occur  to  reset  time  since  disturbance.  First,  if  a  value  less  than  3  years  is 
entered  (or  if  the  age  fields  on  both  the  SHRUBS  and  STDINFO  keywords  are 
left  blank),  initial  time  since  disturbance  will  be  set  to  3  years.  Second,  if  a 
scheduled  thinning  causes  canopy  closure  to  fall  below  50  percent,  and  if  the 
volume  removed  is  20  percent  or  more  of  the  volume  before  thinning,  then  dis- 
turbance time  is  reset  to  3  years.  Finally,  if  the  Regeneration  Establishment 
Model  is  being  used,  and  stand  age  is  reset  to  a  value  less  than  40  years  using 


Calling  the 
Extension 


Overstory  Options 
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Table  2.  — Valid  habitat  type  codes  for  the  SHRUBS  option  (from  Pfister  and  others  1977; 
Steele  and  others  1981) 


Abbreviations 

Codes 

Habitat  types  and  phases 

Pseudotsuga  menziesii  series 

PSME/AGSP 

210 

Pseudotsuga  menziesii/ Agropyron  spicatum 

PSME/FEID 

220 

Pseudotsuga  menziesii/Festuca  idalioensis 

PSME/PHMA 

260 

Pseudotsuga  menziesii/Physocarpus  malvaceus 

PSME/SYAL 

310 

Pseudotsuga  menziesii/Symphoricarpos  albus 

PSME/CARU 

320 

Pseudotsuga  menziesii/Calamagrostis  rubescens 

PSME/CAGE 

330 

Pseudotsuga  menziesii/Carex  geyeri 

PSME/SPBE 

340 

Pseudotsuga  menziesii/Spiraea  betulifolia 

PSME/SYOR 

380 

Pseudotsuga  menziesii/Symptioricarpos  oreophilus 

PSME/ACGL 

390 

Pseudotsuga  menziesii/Acer  glabrum 

PSME/BERE 

395 

Pseudotsuga  menziesii/Berberis  repens 
Abies  grandis  series 

ABGR/SPBE 

505 

Abies  grandis/Spiraea  betulifolia 

ABGR/XETE 

510 

Abies  grandis/Xerophyllum  tenax 

ABGR/COOC 

511 

Abies  grandis/Coptis  occidentalis 

ABGR/VAGL 

515 

Abies  grandis/Vaccinium  globulare 

ABGR/CLUN 

520 

Abies  grandis/Clintonia  uniflora 

ABGR/ACGL 

525 

Abies  grandis/Acer  glabrum 

ABGR/LIBO 

590 

Ahie^  nn^ndi^/l  inn^^p^  tinrp^^li'^ 
Thuja  plicata  series 

THPL/CLUN 

1111     L_l  NwV  ^  Vw/  1  '« 

530 

Thuja  plicata/Clintonia  uniflora 

THPL/ATFI 

540 

Thuja  plicata/ Athyrium  filix-femina 

THPL/OPHO 

550 

Thuja  plicata/Oplopanax  horridum 
Tsuga  heterophylla  series 

TSHE/CLUN 

570 

Tsuga  heterophylla/Clintonia  uniflora 
Abies  lasiocarpa  series 

ABLA/CLUN 

620 

Abies  lasiocarpa/Ciintonia  uniflora 

ABLA/STAM 

635 

Abies  lasiocarpa/Streptopus  amplexifolius 

ABLA/AGGL 

645 

Abies  lasiocarpa/Acer  glabrum 

ABLA/CACA 

650 

Abies  lasiocarpa/Calamagrostis  canadensis 

ABLA/MEFE 

670 

Abies  lasiocarpa/lvlenziesia  ferruginea 

ABLA/XETE 

690 

Abies  lasiocarpa/Xerophyllum  tenax 

ABLA/SPBE 

705 

Abies  lasiocarpa/Spiraea  betulifolia 

TSME/XETE 

710 

Tsuga  mertensiana/Xerophyllum  tenax 

ABLA/VAGL 

720 

Abies  lasiocarpa/Vaccinium  globulare 

ABLA/VAGL 

721 

Abies  lasiocarpa/Vaccinium  globulare,  Vaccinium 
scoparium  phase 

ABLA/VASC 

730 

Abies  lasiocarpa/Vaccinium  scoparium 

ABLA/CAGE 

790 

Abies  lasiocarpa/Carex  geyeri 

ABLA/LUHI 

830 

Abies  lasiocarpa/Luzula  hitchcockii 

the  RE  SET  AGE  keyword,  time  since  disturbance  will  also  be  reset  to  the  same 
value. 

Shrub  statistics  will  be  computed  only  if  a  valid  habitat  type  code  is  encoun- 
tered. Even  if  a  code  other  than  those  listed  in  table  2  has  been  entered  for  the 
stand  on  the  STDINFO  record,  the  user  may  supply  a  separate  habitat  code 
that  affects  only  the  SHRUBS  option.  Thus  shrub  predictions  may  be  made  (at 
the  user's  discretion)  by  substituting  a  similar  valid  habitat  type  code.  Allow- 
ing predictions  to  be  made  for  similar  habitat  types  assumes  that,  even  though 
habitat  type  classification  is  based  on  unique  potential  climax  vegetation,  serai 
community  development  may  not  be  unique  and  similar  habitat  types  may 
respond  with  fairly  similar  shrub  communities  in  the  first  40  years  following 
disturbance.  If  a  substitute  habitat  type  code  is  used,  the  program  writes  a 
cautionary  message  to  that  effect. 
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Shrub  Calibration  If  shrub  information  has  been  recorded  in  the  inventory,  it  may  be  used  to 

Options  adjust  the  embedded  models  to  reflect  unique  variations  in  site  and  environ- 

ment. Shrub  height,  cover,  and  occurrence  models  may  be  cahbrated  using 
either  of  two  types  of  data  collected  according  to  Region  1  Stand  Examination 
Procedures  (USDA  1983).  The  two  methods  are  to  measure  the  height  and  aver- 
age cover  of  up  to  three  distinct  shrub  layers  in  the  stand,  or  alternatively  to 
measure  the  height  and  cover  of  individual  species.  When  real  shrub  measure- 
ments are  provided  to  the  model,  the  information  is  used  to  scale  predictions  to 
match  observed  values.  The  scaling  factors  are  computed  only  once,  at  the 
start  of  the  first  cycle,  and  are  applied  to  the  predictions  for  all  cycles  until  a 
simulated  thinning  occurs.  Because  the  course  of  shrub  development  may  be 
expected  to  be  altered  following  a  thinning,  the  original  information  input  for 
calibration  may  no  longer  be  appropriate.  Thus,  once  a  thinning  occurs,  all  cali- 
bration ceases  and  shrub  predictions  are  no  longer  multiplied  by  the  scaling 
factors. 

SHRBLAYR  Keyword 

This  keyword  is  one  of  two  possible  methods  for  providing  field  data  with 
which  to  calibrate  the  shrub  predictions.  The  SHRBLAYR  keyword  record  con- 
tains six  fields  of  10  columns  each  for  recording  average  height  and  percentage 
of  ground  cover  value  for  each  distinct  shrub  layer  (up  to  three  layers)  in  the 
understory.  There  is  no  inherent  height  or  percentage  of  cover  ranking  of  the 
layers;  they  may  be  entered  in  any  order.  The  information  is  coded  according  to 
the  following  format: 

SHRBLAYR    field  1:  average  height  of  shrub  layer  1,  in  feet 
field  2:  percentage  of  cover  of  shrub  layer  1 
field  3:  average  height  of  shrub  layer  2,  in  feet 
field  4:  percentage  of  cover  of  shrub  layer  2 
field  5:  average  height  of  shrub  layer  3,  in  feet 

field  6:  percentage  of  cover  of  shrub  layer  3 
The  SHRBLAYR  method  of  calibration  sorts  the  shrub  species  by  un- 
calibrated  predicted  height  at  the  beginning  of  the  projection.  Progressing 
down  through  the  species  list  from  predicted  tallest  to  shortest,  the  individual 
uncalibrated  shrub  cover  predictions,  weighted  by  probability  of  occurrence,  are 
summed.  When  the  sum  of  cover  accounts  for  the  same  proportion  as  entered 
for  the  tallest  shrub  layer  on  the  SHRBLAYR  card,  the  summing  ceases,  and 
those  species  are  grouped  into  a  class.  The  process  is  repeated  for  each  input 
layer.  Once  the  classes  are  delineated,  average  predicted  height  and  total  cover 
for  each  layer  are  computed  and  compared  to  the  entered  values  for  the  layers. 
Scaling  factors  are  computed  that  adjust  the  predicted  values  to  match  the 
input  cahbration  values.  The  scaling  factors  are  applied  individually  to  the 
cover  and  height  predictions  of  each  species  within  delineated  classes. 

The  calibration  values  entered  on  the  SHRBLAYR  card  and  the  computed 
scaling  factors  are  displayed  in  the  Shrub  Model  Calibration  Statistics  display 
(fig.  6).  This  display  is  printed  immediately  preceding  the  Shrub  Statistics  dis- 
play if  calibration  is  specified. 
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STAND  ID: 


STAND  GROWTH  PROGNOSIS  SYSTEM  VERSION  5.0  --    INLAND  EMPIRE 

MANAGEMENT  CODE:      NONE  PROGNOSIS  WITH  COVER  EXTENSION  -   SHRUB  LAYER  CALIBRATION 


SHRUB  MODEL  CALIBRATION  STATISTICS 


CALIBRATION  BY  SHRUB  LAYER   ( SHRBLAYR   KEYWORD  CARD) 
AVERAGE  HE IGHT  ( FEET  ) 
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OBSERVED 
VALUES 
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Figure  6.— Shrub  model  calibration  by  shrub  layer,  performed  when  the  following  keyword  records 
are  specified: 

COVER 

SHRUBS  20,0  570.0  3.0 

SHRBLAYR  1.0  20.0  3.0  20.0  6.0  10.0 

END 


SHRUBPC  and  SHRUBHT  Keywords 

The  SHRUBPC  and  SHRUBHT  keywords  are  used  to  supply  calibration  in- 
formation in  instances  where  cover  and/or  average  height  measurements  have 
been  gathered  for  some  or  all  of  the  individual  species. 

The  format  for  the  two  keywords  is  identical,  each  requiring  up  to  five  cards 
for  its  input.  The  leading  card  has  just  the  keyword  entered  on  it.  Following 
are  up  to  four  supplemental  data  records  containing  eight  fields  of  10  columns 
each.  The  first  four  columns  of  each  field  are  used  to  identify  the  shrub  species, 
using  the  abbreviations  given  in  table  1.  The  remaining  six  columns  are  used  to 
enter  the  corresponding  percentage  of  cover  for  the  SHRUBPC  keyword  or 
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height  in  feet  for  the  SHRUBHT  keyword.  For  example,  to  enter  a  50  percent 
cover  statistic  for  Rocky  Mountain  maple,  the  field  would  appear  as: 
ACGL  5.0.^^ 

If  no  data  were  collected  for  a  given  species,  there  should  be  no  entry  for  that 
species.  If  the  species  was  included  in  the  field  survey  but  did  not  occur  in  the 
stand,  it  would  be  entered  in  the  following  manner: 

ACGL  ^^0_ 

Shrub  species  may  be  in  any  order  on  the  supplemental  records.  Enter 
"  —  999"  in  the  shrub  code  field  following  the  last  shrub  entry  to  signify  the  end 
of  the  SHRUBPC  (or  SHRUBHT)  data. 


STAND  GROWTH   PROGNOSIS  SYSTEM  VERSION  5.0  --    INLAND  EMPIRE 

STAND   ID:      521)8112  MANAGEMENT  CODE:      NONE  PROGNOSIS  WITH  COVER  EXTENSION  -  SHRUB  SPECIES  CALIBRATION 


SHRUB  MODEL  CALIBRATION  STATISTICS 


CALIBRATION  BY    INDIVIDUAL   SPECIES   (SHRUBHT  AND/OR  SHRUBPC   KEYWORD  CARDS): 

SHRUB  HEIGHT  (FEET)  PERCENT  COVER 


SHRUB 

OBSERVED 

PREDICTED 

SCAL 1 NG 
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PREDICTED 

SCALI NG 

SPEC  1 ES 
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Figure  7.— Shrub  model  calibration  by  species,  performed  when  the  following  keyword  records 
are  specified: 


COVER 
SHRUBS 
SHRUBHT 
ACGL  10 
ARUV  0 
SHRUBPC 
ACGL  2 
ARUV  0 
END 


OME  F  E 
0  VASC 


OMEFE 
0  VASC 


2  0.0 


0  VAME 
OCARX 

0  VAME 
OCARX 


5  7  0.0 


0  PAMY 
0-999 


0  P  AMY 
0-999 


3  .  0 

2.0LIBO  0.5PREM  O.OXETE  O.OFERN  2.0 
6.0LIBO       10    OPREM         O.OXETE         O.OFERN  5.0 
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The  SHRUBPC  and  SHRUBHT  keywords  do  not  both  have  to  be  present  in 
a  given  projection,  although  it  is  desirable.  If  only  one  of  the  keywords  is  pres- 
ent, only  that  portion  of  the  calibration  will  be  performed.  The  SHRBLAYR 
keyword  should  not  be  included  when  using  the  SHRUBPC  and/or  SHRUBHT 
keywords. 

The  data  supplied  on  the  SHRUBPC  and  SHRUBHT  keyword  cards  are  used 
to  adjust  probability  of  occurrence,  height,  and  cover  predictions  for  individual 
shrub  species.  Scaling  factors  are  computed  for  each  species  as  the  ratio  of  ac- 
tual to  predicted  height  and  cover  at  the  start  of  the  projection.  Any  species 
recorded  as  absent  is  given  a  zero  probability  of  occurrence.  In  each  cycle, 
these  scaling  factors  are  applied  to  the  appropriate  species  and  prediction.  Scal- 
ing factors  are  ignored  after  the  first  simulated  thinning. 

Height  and  cover  calibration  factors  by  species  are  output  in  the  Shrub 
Model  Calibration  Statistics  Display  (fig.  7). 

Additional  There  are  seven  additional  keywords  in  the  COVER  extension.  The  first. 

Keywords  SHOWSHRB,  is  used  to  select  up  to  six  understory  species  for  which  output 

will  always  be  displayed.  These  are  in  addition  to  the  nine  species  that  account 
for  the  most  cover.  It  requires  one  supplemental  data  record  containing  shrub 
species  codes  (table  1)  in  six  fields  of  10  columns  each.  The  four-character  codes 
must  be  right-justified  within  the  fields. 
Card  1:  SHOWSHRB 

Card  2:  Cols.  1  —  10:  abbreviation  for  first  species 


Cols.  51  —  60:  abbreviation  for  last  species 
Three  keywords  are  used  to  turn  off  printing  of  the  displays: 
NOCOVOUT    Suppress  output  of  Canopy  Cover  Statistics  display; 

default  =  table  printed. 
NOSHBOUT    Suppress  output  of  the  Shrub  Statistics  display; 

default  =  table  printed. 
NOSUMOUT   Suppress  output  of  the  Canopy  and  Shrubs  Summary  display; 
default  =  table  printed. 
The  final  three  keywords  are  used  to  request  information  about  the  program, 
print  intermediate  debug  information,  and  to  signify  the  end  of  the  COVER 
options: 

DATELIST   Print  date  of  last  revision  of  COVER  model  subprograms 

and  common  areas. 
DEBUG        Request  printout  of  the  results  of  most  calculations  for 

all  tree  and  shrub  records  (caution:  voluminous  output!). 

field  1:  Cycle  number  in  which  debug  output  is  to  be  printed; 
default  =  output  printed  in  all  cycles. 
END  Signify  the  end  of  COVER  keywords  and  return  control 

to  the  main  program. 

The  COVER  keywords  for  the  example  stand  are  echoed  in  the  Prognosis 
Model  Input  Summary  table  (fig.  8). 
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CANOPY 
Submodels 


OVERVIEW  OF  THE  COVER  SUBMODELS 

Crown  Width  and  Stand  Canopy  Closure.— COVER  predicts  crown  develop- 
ment for  the  11  conifer  species  listed  in  table  3.  Logarithmic  regression  equa- 
tions are  used  to  predict  individual  tree  crown  width  from  species,  d.b.h., 
height,  and  crown  length  for  trees  3.5  inches  d.b.h.  and  larger.  For  trees  less 
than  3.5  inches  d.b.h.,  crown  width  is  a  function  of  species,  height,  crown 
length,  and  stand  basal  area  (Moeur  1981).  Coefficient  values  for  the  crown 
width  models  for  large  and  small  trees  are  shown  in  tables  4  and  5, 
respectively. 

Individual  tree  crown  area  is  computed  as  the  area  of  a  circle  with  diameter 
equal  to  predicted  crown  width.  Stand  canopy  closure  is  computed  from  the 
sum  of  the  tree  crown  areas, 

L  crown  areas  (ft^/acre) 


canopy  closure  = 


43,560  ft^/acre 


X  100  percent. 


Table  3. 

—Tree  species  recognized  by  COVER 

Code 

Common  name 

Scientific  name 

WP 

Western  wliite  pine 

Pinus  monticola 

L 

Western  larch 

Larix  occidentalis 

DF 

Douglas-fir 

Pseudotsuga  menziesii 

GF 

Grand  fir 

Abies  grandis_ 

WH 

Western  hemlock 

Tsuga  heterophylla 

C 

Western  redcedar 

Thuja  pi  i  cat  a 

LP 

Lodgepole  pine 

Pinus  contorta 

S 

Engelmann  spruce 

Picea  engelmannii 

AF 

Subalpine  fir 

Abies  lasiocarpa 

PP 

Ponderosa  pine 

Pinus  ponderosa 

Other 

Whitebark  pine 

Pinus  albicaulis 

Table  4.— Coefficients  for  estimating  crown  width  of  trees 
3.5  inches  d.b.h.  and  larger  (Moeur  1981): 
ln(crown  width)  =  b^  +  bJn{D)  +  bjInlH)  +  bjInlCL) 


Variable  coefficients^ 


Intercept 

ln(H) 

Species^ 

bo 

b2 

WP 

4.30800 

-1.37264 

L 

2.31359 

-  .80919 

DF 

3.02271 

-  1.00486 

GF 

2.20611 

-  .76936 

WH 

1.32772 

-  .52554 

C 

2.79784 

-  .89666 

LP 

1.06804 

-  .55987 

S 

3.76535 

-1.18257 

AF 

1.74558 

-  .73972 

PP 

1.62365 

-  .68098 

Other 

-  .91984 

-  .07299 

Variables       Variable  coefficients 

ln{D)  bi  =  1.08137 

In(CL)  b3  =  .29786 


'Species  codes  are  given  in  table  3. 
^Definition  of  variables: 

D  =  diameter  breast  height  (inches) 

H  =  tree  height  (ft) 

CL  =  crown  length  (ft). 


21 


Table  5.— Coefficients  for  estimating  crown  width  of  trees 
less  than  3.5  inches  (Moeur  1981): 
ln(crown  width)  =  b,ln(H)  +  b2ln(CL)  +  bglniBA) 

Variable  coefficients^ 


ln(H) 


Species^ 

WP 

0.37031 

1 

DF 

.32874 

GF 

.38503 

WH 

.25622 

C 

.46452 

LP 

.26342 

S 

.33089 

AF 

.33722 

PP 

.36380 

Other 

.07049 

Variables 

Variable  coefficients 

In(CL) 
In(BA) 

b2  =  0.28283 
b3  =  .04032 

''Species  codes  are  given  in  table  3. 
^Definition  of  variables: 

H  =  tree  hieight  (ft) 

CL  =  crown  length  (ft) 

BA  =  stand  basal  area  (ft^/acre). 


In  the  following  discussion,  model  behavior  is  displayed  in  the  plots  of  simu- 
lation results  from  five  stands  whose  site  characteristics  are  Usted  in  table  6. 
In  each  100-year  simulation,  the  stand  was  clearcut  in  period  1,  the  site  was 
prepared  by  burning,  and  a  new  tree  list  predicted  using  the  Regeneration 
Establishment  Model  (Ferguson  and  Crookston  1984).  These  "bare-ground" 
regeneration  projections  were  used  to  compare  responses  of  the  crown  relation- 
ships to  changes  in  stand  structure  and  density  through  time. 

Table  6.— Site  characteristics  of  the  stands  used  to  examine  crown  model 
behavior.  In  each  simulation,  the  stand  was  clearcut,  the  site  was 
prepared  by  burning,  and  regenerated  using  the  Regeneration 
Establishment  Model 

Code      Stand       Location        Habitat  type      Aspect     Slope  Elevation 


Percent  Feet 


E 

S248112 

St.  Joe 

570  (TSHE/CLUN) 

NW 

30 

3,400 

W 

Weippe 

Clearwater 

530  (THPL/CLUN) 

NW 

20 

4,000 

C 

Cranberry 

Clearwater 

530  (THPL/CLUN) 

S 

10 

3,000 

G 

Grouse 

Clearwater 

520  (ABGR/CLUN) 

N 

10 

3,100 

S 

Silver 

Clearwater 

520  (ABGR/CLUN) 

S 

20 

3,000 

Canopy  closure  follows  an  increasing  sigmoidal  pattern  over  time  on  regener- 
ated stands  (fig.  9).  Cover  increases  fairly  rapidly  between  0  and  20  years,  as 
the  first  and  second  waves  of  regeneration  produced  by  the  establishment 
model  enter  the  tree  list.  Then,  even  though  numbers  of  trees  start  to  decline 
beyond  20  years,  canopy  closure  increases  as  individual  tree  crown  width  con- 
tinues to  expand.  In  the  stands  where  canopy  closure  reaches  100  percent  and 
greater,  crown  cover  peaks  and  then  dechnes  beyond  about  50  years  in  the 
projection  (stands  E,  W,  and  G).  Canopy  closure  is  incomplete  on  the  stands 
where  establishment  is  poor  (C  and  S).  Here,  crown  cover  levels  off  about  70 
years  after  regeneration,  rather  than  peaking  and  declining. 
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BARE-GROUND  PROJEC^ONS  ON  USER'S  MANUAL  EXAMPLE  STANDS.  TABLE  6 


(a) 


LEGENDi  STRND 


TIKE  SINCE  nECENEnmiON  IN  TERB3 


 GROUSE 


Figure  9.— Crown  cover  predictions  for  the  regener- 
ated stands  in  table  6:  (a)  predicted  trees  per  acre 
versus  time;  (b)  canopy  ciosure  versus  time;  (c)  aver- 
age tree  crown  widtli  versus  time. 
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In  stands  where  initial  densities  exceed  about  400  trees  per  acre  (E.  W,  G), 
average  crown  width  increases  continually  through  time  as  the  number  of  trees 
in  the  stand  decreases.  Presumably,  this  is  related  to  the  amount  of  mortality 
and  to  the  types  of  trees  that  are  dying.  As  number  of  trees  decreases,  crown 
width  on  surviving  trees  continues  to  increase  through  time.  In  addition,  more 
of  the  mortality  in  the  older  stands  is  accounted  for  by  suppressed  trees  in  the 
understory,  leaving  trees  with  relatively  more  vigorous  crowns.  Below  initial 
densities  of  400  trees  per  acre,  there  is  little  relationship  between  individual 
tree  crown  width  and  number  of  trees,  or  between  percentage  of  cover  and 
trees.  Crown  width  increases  throughout  the  length  of  the  projection  in  the 
understocked  stands,  even  though  tree  numbers  remain  relatively  constant,  in- 
dicating that  stand  density  is  not  great  enough  to  be  limiting  to  individual 
crown  development.  Under  apparent  lack  of  competition  for  space,  the  tree 
crowns  in  stands  C  and  S  are  free  to  expand.  In  general,  the  model  predicts 
that  crowns  will  be  wider  in  stands  that  start  out  with  more  open  conditions. 

Figure  10  shows  the  relationship  of  overstory  cover  to  stand  basal  area  and 
the  effect  of  competition  between  crowns.  Canopy  closure  is  monotonically  in- 
creasing with  increasing  basal  area  in  all  stands.  Canopy  closure  attains  a  max- 
imum and  then  declines  in  the  fully  stocked  stands,  which  reach  maximum 
basal  area  and  then  level  off  (E,  W,  and  G).  Canopy  closure  levels  off  but  does 
not  begin  to  decline  on  stands  that  do  not  reach  maximum  basal  area  within 
the  time  frame  of  the  projection  (C  and  S).  For  a  given  basal  area  level,  crown 
width  is  greater  in  the  more  open  stands,  this  difference  increasing  throughout 
the  projection. 

BARE-GROUND  PROJECTIONS  ON  USER'S  MANUAL  EXAMPLE  STANDS,  TABLE  6 


350- 


Figure  10.— (a)  Stand  basal  area  over  time;  (b) 
predicted  canopy  closure  versus  basal  area:  (c)  aver- 
age tree  crown  width  versus  basal  area  for  the 
regenerated  stands  in  table  6. 
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(C)  BftSRL  RREH  IN  Sa.n.  PER  ACRE 


Fig  10.— (Con.) 
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Crown  Foliage  Biomass.— Foliage  biomass  is  predicted  for  individual  trees 
using  logarithmic  linear  regression  equations  (Moeur  1981).  As  in  the  crown 
width  functions,  there  are  different  equation  forms  for  trees  less  than  3.5  inches 
d.b.h  and  those  3.5  inches  d.b.h.  and  larger.  The  equations  predict  foliage  bio- 
mass from  species,  d.b.h.,  change  in  squared  diameter,  height,  crown  length, 
trees  per  acre,  and  relative  diameter  (d.b.h. /quadratic  mean  stand  diameter)  for 
large  trees,  and  from  species,  height,  crown  length,  and  trees  per  acre  for  small 
trees.  Coefficients  for  the  models  are  Usted  in  tables  7  and  8. 

Table  7.— Coefficients  for  estimating  foliage  biomass  of  trees 
3.5  inchies  and  larger  (fVloeur  1981):  In(biomass)  = 
bg  +  b^ln(D)  +  bjIniH)  +  b3ln(CL)  +  bJniDDS)  + 
b5ln{TPA)  +  b5ln(DREL) 


Variable  coefficients^ 

Intercept 

Species^ 

bo 

WP 

2.66607 

L 

1.75654 

DF 

2.70587 

GF 

3.11508 

WH 

2.65457 

C 

3.05935 

s 

3.30085 

AF 

3.06017 

PP 

2.45249 

Otfier 

2.62251 

Variables 

Variable  coefficients 

ln(D) 

=  1.468547 

ln(H) 

=  1.07705 

In(CL) 

=  .69082 

In(DDS) 

b^  =  .30885 

In(TPA) 

bg  =  .14210 

In(DREL) 

bg  =  .39924 

''Species  codes  are  given  in  table  3. 
^Definition  of  variables: 

D  =  diameter  breast  height  (inches) 

H  =  tree  height  (ft) 

CL  =  crown  length  (ft) 

DDS  =  change  in  squared  diameter  (in^) 

TPA  =  trees  per  acre 

DREL  =  d.b.h/quadratic  mean  d.b.h. 
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Table  8.— Coefficients  for  estimating  foliage  biomass  of  trees 
less  than  3.5  incties  d.b.h.  (Moeur  1981): 
In(biomass)  =  bg  +  bJn(H)  +  b2ln(CL)  +  b3ln(TPA) 


Variable  coefficients^ 


Intercept 

In(CL) 

Species^ 

bo 

b2 

WP 

- 1.94951 

1.22023 

L 

-4.73762 

1.98479 

DF 

-2.05828 

1.25837 

c 
or 

-  2.43200 

1 .60270 

vv  rn 

A  17yic;c 
—  ^.  \  (  HoKy 

c 

-2.24876 

1.37600 

LP 

-3.13488 

1.62368 

S 

-2.93508 

1.96125 

AF 

-1.60998 

1.32649 

PP 

-2.74410 

1.58171 

Otfier 

-1.63387 

1.35092 

Variables 

Variable  coefficients 

ln(H) 

b^  =  0.40350 

In(TPA) 

bj  =  .12975 

^Species  codes  are  given  in  table  3. 
^Definition  of  variables: 

H  =  tree  height  (ft) 

CL  =  crown  length  (ft) 

TPA  =  trees  per  acre. 


In  general,  the  response  patterns  of  total  stand  foliage  biomass  to  trees  per 
acre  and  basal  area  over  time  are  similar  to  those  for  canopy  closure  (fig.  11). 
Total  foliage  is  greater  in  the  stands  starting  out  with  higher  initial  densities, 
peaking  at  50  to  60  years,  and  then  declining  (stands  E,  W,  and  G).  Stand  foli- 
age does  not  decline  as  sharply  as  does  percentage  of  canopy  closure  near  the 
end  of  the  projection.  Stands  with  low  initial  densities  show  gradually  increas- 
ing values  of  foliage  throughout  the  projection  (C  and  S).  Foliage  production 
also  depends  on  the  species  composition  of  the  mature  stands.  Stands  W  and  C 
have  higher  proportions  of  mature  trees  in  grand  fir,  cedar,  and  Douglas-fir. 
These  three  species  have  greater  predicted  foliage  values  for  a  given  set  of 
stand  conditions.  Individual  tree  foliage  development  patterns  through  time  are 
quite  similar  to  those  discussed  for  crown  width. 
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BARE-GROUND  PROJECTIONS  ON  USER'S  MANUAL  EXAMPLE  STANDS.  TABLE  6 
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Figure  1 1.— (Con. 
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Crown  Shape.— Individual  tree  crown  shape  is  predicted  each  projection  cycle 
using  a  linear  discriminant  function  (Moeur  1983).  Tree  crowns  are  classified 
into  one  of  five  shapes— circular,  triangular,  neiloid,  parabolic,  or  elliptic— using 
species,  height,  d.b.h.,  crown  length,  crown  radius,  crown  ratio,  and  trees  per 
acre  as  discriminating  variables.  All  tree  crowns  are  assumed  to  have  a  circular 
bottom  (fig.l2). 

Crown  shape  is  used  in  three  places  in  COVER.  First,  crown  volume  within 
vertical  height  classes  in  the  stand  is  computed  by  summing  sections  of  in- 
dividual crowns,  using  formulas  to  integrate  the  five  different  sohds  of  revolu- 
tion in  figure  12.  Second,  crown  profile  area  within  height  classes  is  computed 
by  summing  the  lateral  area  of  individual  crown  sections  (fig.  3b).  Finally,  tree 
foliage  biomass  is  distributed  within  height  classes  by  the  proportion  of 
frustum  volume  within  height  classes  to  total  crown  volume.  An  inner  sene- 
scent cone  is  not  considered;  that  is,  foliage  is  assumed  to  be  distributed  uni- 
formly throughout  the  crown. 


PROFILE  SHAPE 


PROFILE  SHAPE 

PARABOLA  ELLIPSE 


SOLID  FORM  PARABOLOID  ELLIPSOID 

Figure  12.— Individual  crowns  are  classified  into  one  of  five  shapes  in  the  CANOPY 
option  (after  Mawson  and  others  1976). 
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The  equations  for  predicting  understory  species  relationships  have  been 
synthesized  from  the  work  of  different  people  at  different  times.  An  early 
version  of  the  SHRUBS  extension,  known  as  BROWSE  (Scharosch  1980), 
incorporated  studies  on  shrub  development  in  grand  fir,  cedar,  and  hemlock 
habitat  types  conducted  by  Irwin  and  Peek  (1979).  In  Prognosis  5.0,  these 
relationships  have  been  replaced  by  probability  of  occurrence  equations  for 
individual  understory  species  developed  by  Scharosch  (1984),  and  height  and 
cover  equations  developed  by  Laursen  (1984).  Irwin  and  Peek's  work  was  based 
on  a  subset  (about  2,200  plots)  of  the  data  used  by  Laursen  (1984)  and 
Scharosch  (1984).  The  expanded  data  (about  10,000  plots  located  in  about  500 
stands)  include  measurements  from  Douglas-fir  and  subalpine  fir  ecosystems  as 
well  as  the  original  grand  fir-cedar-hemlock  types.  Table  1  lists  the  species, 
range  of  habitat  types,  and  sources  of  information  for  which  predictions  are 
currently  made. 

Inside  the  COVER  program,  the  probability  of  any  shrub  cover  on  the  site, 
and  total  shrub  cover  given  that  the  probability  is  greater  than  zero,  are 
predicted  first.  Then,  probability  of  occurrence  is  calculated  for  each  species 
individually  using  total  shrub  cover  as  a  predictor.  Next,  heights  are  predicted 
for  each  species,  also  using  total  shrub  cover.  The  species  are  sorted  in  order 
from  tallest  to  shortest  predicted  height.  Then,  progressing  down  through  the 
sorted  list,  individual  species  cover  is  calculated  using  predicted  species  height 
and  the  amount  of  overtopping  cover.  Species  cover  is  weighted  by  species 
probability  of  occurrence.  Finally,  the  cover  values  are  summed  and  reported  as 
total  understory  cover  for  the  plot. 

Total  Shrub  Cover.— Predictions  for  probability  of  any  shrubs,  and  for  total 
shrub  cover,  conditional  on  probability,  are  taken  from  Laursen  (1984).  The 
probability  that  shrub  cover  exists  given  the  described  stand  conditions  is 
calculated  and  reported  in  the  summary  display.  This  value  expresses  the 
proportion  of  1/300-acre  plots  on  which  shrub  cover  is  expected  to  be  greater 
than  zero.  It  is  computed  from  a  logistic  regression  model  using  slope, 
elevation,  overstory  basal  area,  habitat  type,  time  since  disturbance,  and  the 
interaction  of  time  since  disturbance  and  basal  area.  Next,  total  shrub  cover  is 
predicted  using  a  lognormal  linear  regression  model  fit  on  plots  in  the  original 
data  where  cover  was  actually  present.  Additional  variables  in  this  model  are 
type  of  disturbance  and  the  time  and  type  of  disturbance  interaction.  Total 
shrub  cover  predicted  by  Laursen's  model  for  three  hypothetical  treatments  is 
shown  in  figure  13. 
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cavEB 


TOTAl.  PERCENTAGE  SHRUB  COVER 
DISTURBflNCE  TTPE=NONE 


ISO  - 


100 


Figure  13.— Predicted  total  percent  stirub 
cover  relative  to  time  since  disturbance  and 
overstory  basal  area  following  (a)  no  site 
preparation,  (b)  mectianical  disturbance,  and 
(c)  burning.  Variables  tield  constant  are  slope 
=  0.25,  elevation  =3.500  ft,  fiabitat  type  = 
ABGR/CLUN  (from  Laursen  1984). 
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Table  9.— Variables  used  to  predict  probability  of  occurrence,  height,  and  per- 
centage of  cover  for  the  shrub  species  listed  in  table  1 

Continuous  variables 


Overstory  basal  area  (ft^/acre) 

Stand  elevation  (100's  of  feet) 

Elevation  squared  (10,000's  of  feet) 

Total  percent  shrub  cover 

Slope  (percent/100) 

Slope  X  sin(aspect) 

Slope  X  cos(aspect) 

Time  since  site  disturbance  (years) 

Overstory  basal  area  x  time  since  disturbance 


Categorical  variables 

Overstory  climax  species  (includes  habitat  type):^ 

Douglas-fir  (210,  220,  260,  310,  320,  330,  340,  380,  390,  395) 
Grand  fir  (505,  510,  511,  515,  520,  525,  590) 
Western  redcedar  (530,  540,  550) 
Western  hemloci^  (570) 

Subalpine  fir/mountain  hemlock  (620,  635,  645,  650,  670,  690,  705,  710, 

720,  721,  730,  790,  830) 
Understory  climax  union  (includes  habitat  type):'' 

ABGR/CLUN,  COOC,  XETE,  LIBO,  ABLA/STAM,  ABLA/LUHI  (510,  511,  520,  590, 

635,  690,  710,  830) 
THPL  series  (530,  540,  550) 

Tall  shrub  (260,  390,  515,  525,  645,  670,  720,  721,  730) 
Low  shrub  (310,  340,  380,  395,  505,  705) 
TSHE/CLUN,  ABLA/CLUN  (570,  620) 
Grasses  (210,  220,  320,  330,  650,  790) 
National  Forest  grouping: 
Boise,  Payette 
Nezperce 

Clearwater,  Coeur  d'Alene,  Lolo,  St.  Joe 
Panhandle,  Colville,  Kaniksu,  Kootenai 
Physiography: 
Bottom 
Lower  slope 
Midslope 
Upper  slope 
Ridge 

Type  of  site  disturbance: 
None 

Mechanical 

Burn 

Road 

Type  of  disturbance  x  overstory  basal  area 
Type  of  disturbance  x  time  since  disturbance 


'See  table  2  for  habitat  type  code  definitions. 


Species  Probability  of  Occurrence.— Predicted  probability  of  occurrence 
equations  for  the  species  in  table  1  are  from  Scharosch  (1984).  He  uses  a 
logistic  multiple  regression  model  to  predict  species  occurrence  from  the 
continuous  and  categorical  variables  in  table  9.  The  logistic  model  produces  a 
sigmoidal  curve  with  predicted  values  restricted  to  the  closed  interval  (0,1). 
Representative  responses  of  predicted  probability  of  occurrence  are  plotted  in 
figure  14. 
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CEVE  -  PROBABILITY  VS.  BASAL  AREA.  BY  TIME  SINCE  DISTURBANCE 

VflBIflBLES  HELD  CONSTflNTi  SHRUB  CI3VER=5(r/:.  Q.EV=iJ0OO,  RSPECT=180, 
SLflP£=.30.  HfiB=GF/CLUN.  PKTS=BOnOM.  OIST=MECHflNICflL.  NF=.ST.JOE 
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PHMA  -  PROBABILfTY  VS.  SLOPE  AND  ASPECT,  BY  HABITAT  TYPE 
VHRIflBLES  HELD  CONSTflNTi  SHRUB  C0VER=50X,  ELEV=14000, 
BH=20.  TIME=20.  PHT3=MI0SL0PE.  0I3T=NBNE.  NF=ST.JOE 


SLOPE 

(b) 

Figure  14.— (a)  Predicted  probability  of  occurrence  across  a  range 
of  overstory  basal  area  and  time  since  disturbance  for  Ceanothus 
velutinus;  (b)  predicted  probability  of  occurrence  by  slope  and 
aspect  for  Physocarpus  malvaceus  on  PSI^E/PHI^A  and 
ABGR/CLUN  habitat  types:  (c)  predicted  probability  of  occurrence 
by  disturbance  type  for  Ceanothus  sanguineus  and  Acer  glabrum 
(from  Scfiarosch  1984). 
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ACGL  AND  CESA  -  PROBABILITY  VS.  TIME,  BY  DISTURBANCE  TYPE 

VflHIFIBLES  HELD  CONSTRNTt  SHRUB  OmB^SfT/..  ELEV=1JOOO,  Bfi=20, 
SL0PE=.20,  fiSPECT=22S,  HflB=HRC/CLUN,  PHYS=RIDGE.  NF=ST.jaE 


(c) 


LEGEND!  DISTURB 


15  20  25 

TIME  SINCE  DISTURBfiNCE  (TEARS) 


—  RCGL  -  BURN 

—  CESfl  -  BURN 


RCGL  -  MECH 
CESR  -  MECH 
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RCGL  -  NONE 
CESfl  -  NONE 


Fig  i4.  — (Con. 


Individual  Species  Height.— Laursen  (1984)  uses  either  Hnear  or  lognormal 
regression  model  forms  to  predict  the  average  heights  of  individual  species  in 
the  described  stand.  Model  forms  differ  between  species,  but  most  contain  time 
since  disturbance,  overstory  basal  area,  and  predicted  total  shrub  cover  and  its 
residual  (the  difference  between  observed  and  predicted  values  when  an  ob- 
served value  is  supplied  by  the  user)  as  independent  variables.  Various  transfor- 
mations of  the  independent  variables  in  table  9  are  included  to  represent  other 
site  and  treatment  effects.  The  equations  are  detailed  in  Laursen 's  paper. 
Representative  response  patterns  of  predicted  height  over  time  and  stand  basal 
area  are  plotted  in  figure  15. 

Individual  Species  Cover.— Percentage  of  cover  by  species  follows  lognormal 
or  logistic  distributions,  conditional  on  the  presence  of  the  species  in  the  stand 
(Laursen  1984).  Cover  for  most  species  is  a  function  of  predicted  species  height 
and  its  residual  (observed  minus  predicted  when  observed  values  are  supplied), 
overtopping  by  taller  species  (the  percentage  of  predicted  cover  above  current 
height),  time  since  disturbance,  type  of  disturbance,  overstory  basal  area,  and 
site  conditions.  Representative  plots  of  the  response  of  species  cover  to  basal 
area  over  time  are  shown  in  figure  16. 

Twig  Production  and  Dormant  Season  Shrub  Biomass.— Total  current  year's 
twig  production  in  twigs  per  square  foot  and  total  dormant  season 
aboveground  shrub  biomass  in  pounds  per  acre  are  computed  only  for  the  16 
species  noted  in  table  1  and  only  on  ABGR/CLUN,  THPL/CLUN,  and 
TSHE/CLUN  habitat  types  (Irwin  and  Peek  1979).  Twig  production  is  a  log- 
linear  regression  equation  dependent  on  time  since  disturbance,  overstory 
crown  competition  factor  (CCF)  (Krajicek  and  others  1961),  and  habitat  type. 
Shrub  biomass  is  also  log-linear,  predicted  from  time  since  disturbance  and  CCF. 
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Figure  15.— Predicted  shrub  height  relative  to  overstory  basal 
area  and  time  since  disturbance  for  (a)  Acer  glabrum,  (b) 
Ceanothus  sanguineus,  and  (c)  Vaccinium  scoparium.  Variables 
held  constant  are  slope  =0.25.  aspect  =  east,  elevation  =  3.500 
ft.  habitat  type  =  ABGR/CLUN.  physiography  =  midslope.  distur- 
bance type  =  burn,  forest  =St.  Joe  (from  Laursen  1984). 
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Figure  16.— Predicted  shrub  cover  relative  to  overstory  basal  area 
and  time  since  disturbance  for  (a)  Acer  glabrum.  (b)  Ceanothus 
sanguineus,  and  (c)  Vaccinium  scoparium.  Variables  held  constant 
are  slope  =  0.25.  aspect  =  east,  elevation  =  3.500  ft.  habitat  type 
=  ABGR/CLUN.  physiography  =  midslope.  disturbance  type  =  burn, 
forest  =  St.  Joe  (from  Laursen  1984).  Values  shown  are  for  predicted 
cover  prior  to  multiplication  by  probability  of  occurrence. 
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Stand  Successional  One  subprogram  within  the  COVER  extension  computes  and  displays  a 

Stage  "stand  successional  code"  that  is  related  to  the  vertical  structure  of  both 

shrubs  and  trees  (Peterson  1982).  The  purpose  of  the  classification  is  to  provide 
a  basis  for  relating  wildlife  use  to  a  particular  type  of  stand  that  now  exists,  or 
that  will  result  from  management.  The  codes  are  listed  in  table  10. 

The  classification  is  a  function  of  time  since  stand  disturbance,  crown  compe- 
tition factor  (CCF),  average  tree  height,  a  selectively  defined  "average"  stand 
diameter,  and  average  tall  shrub  height  (average  predicted  height,  weighted  by 
predicted  cover  for  Acer  glabrum,  Alnus  sinuata,  Amelanchier  alnifolia, 
Ceanothus  sanguineus,  Ceanothus  velutinus,  Cornus  stolonifera,  Holodiscus  dis- 
color, Prunus  emarginata,  Prunus  virginiana,  Salix  spp.,  Sambucus  spp.,  and 
Sorbus  spp.). 

To  compute  average  stand  diameter,  a  series  of  logical  tests  determines 
whether  the  stand  is  even-aged,  two-storied,  or  all-aged  based  on  the  distribu- 
tion of  trees  per  acre  and  percentage  of  cover  by  10-ft  height  classes.  For  even- 
aged  stands  (a  stand  in  which  90  percent  of  the  total  canopy  closure  is 
accounted  for  by  trees  within  a  30-ft  height  range),  the  root  mean  square  di- 
ameter of  the  stand  is  used  as  the  average  diameter.  For  two-storied  stands 
(the  two  most  dense  20-ft  layers  must  be  separated  by  20  ft  or  more),  the  root 
mean  square  diameter  for  the  most  dense  20-ft  layer  only  is  used  as  the  aver- 
age diameter.  For  all-aged  stands  (the  most  dense  10-ft  layer  contains  less  than 
20  percent  of  the  total  canopy  closure,  the  three  most  dense  10-ft  layers  contain 
less  than  50  percent,  etc.),  the  average  diameter  is  taken  to  be  the  root  mean 
square  diameter  of  the  three  most  dense  10-ft  layers.  The  stand  is  then  classi- 
fied according  to  the  scheme  in  table  10.  Note  that  restrictions  for  stages  1  to 
4  are  of  the  type  "CCF  less  than  30  or  average  tree  height  less  than  1  ft."  For 
stages  5  and  6,  there  is  a  "percent  shrub  cover"  or  "average  tall  shrub  height" 
restriction. 


Table  10.— Classification  scheme  for  assigning  stand  successiona 

stage  code  (after  Peterson  1982) 

Tall  shrub 

Tall  shrub 

Tall  shrub 

Recent 

Low 

Medium 

with  no 

with  few 

with  mostly 

Sapling 

Pole 

Mature 

Old-growth 

disturbance 

shrub 

shrub 

conifers 

conifers 

conifers 

timber 

timber 

timber 

timber 

Condition 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

Time  since  stand 

disturbance  (years) 

<5 

Crown  competition 

<30 

<30 

<30 

<30 

30-50 

50-100 

>100 

>100 

>125 

>100 

factor 

or 

or 

or 

or 

Average  tree  hieigtnt 

<1 

<1 

<1 

<1 

<5 

>5 

(ft) 

Percent  slirub  cover 

<25 

>70 

>50 

>30 

or 

or 

Average  "tall  shrub" 

height  (ft) 

<1.0 

<2.5 

2.5-5.0 

>5.0 

>5.0 

>5.0 

"Average"  stand 

diameter  (inches) 

<4 

4-11 

11-24 

>24 
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USING  THE  COVER  EXTENSION  AS  A  MANAGEMENT 
TOOL 

The  general  comments  in  this  section  are  intended  to  guide  the  user  in 
applying  information  produced  by  the  combined  Prognosis/COVER  model  in  a 
broader  planning  context.  There  is  little  specific  COVER  output  that  directly 
interprets  wildlife,  hydrologic,  or  insect  pest  relationships.  Instead,  it  is  up  to 
the  user  to  interpret  the  information  specific  to  his  or  her  application.  An 
important  criterion  in  developing  COVER  was  to  make  it  broad  enough  in 
design  for  many  applications,  but  primarily  to  link  vegetation  changes  with 
nontimber  resources.  COVER  can  be  a  useful  tool  for  decisionmaking  when 
combined  with  knowledge  of  a  specific  resource  ecology  and  its  relation  to 
vegetation  management  systems. 

Wildlife  Habitat  Many  of  the  shrub  and  tree  cover  development  values  produced  by  COVER 

Applications  can  be  related  to  wildlife  habitat.  An  example  stand  projection  illustrates  how 

displays  generated  from  COVER  values  can  be  used  to  compare  vegetation 
changes  and  alternative  treatment  effects  on  wildlife  habitat.  The  example 
presented  is  a  stand  that  was  inventoried  in  1984  at  145  years  of  age.  Initial 
stand  density  is  459  trees  per  acre,  composed  of  an  understory  of  Engelmann 
spruce  beneath  a  sparse  overstory  of  Douglas-fir,  western  larch,  and  grand  fir. 
The  stand  is  on  the  Nezperce  National  Forest,  ABGR/VAGL  habitat  type, 
northeast  aspect,  50  percent  slope,  at  5,800  ft  elevation.  Two  silvicultural  treat- 
ments were  simulated  and  compared.  In  the  first,  the  lower  and  upper  portions 
of  the  diameter  distribution  were  removed  in  cycle  1,  leaving  60  Douglas-fir, 
larch,  and  spruce  trees  per  acre,  with  diameters  between  18  and  25  inches.  At 
the  end  of  the  second  cycle,  natural  regeneration  was  predicted  to  be  780  trees 
per  acre,  composed  of  60  percent  grand  fir,  25  percent  Douglas-fir,  and  15  per- 
cent spruce.  Twenty  years  after  the  initial  thin,  all  trees  greater  than  7  inches 
d.b.h.  were  removed,  and  the  regenerated  stand  was  grown  to  age  100.  This 
prescription  is  referred  to  as  "two-step  shelterwood. "  The  second  prescription, 
"clearcut,"  cut  all  trees  in  the  stand  in  cycle  1.  Natural  regeneration  at  the  end 
of  the  second  10-year  cycle  was  450  trees  per  acre,  of  which  65  percent  were 
grand  fir,  30  percent  Douglas-fir,  and  5  percent  larch.  This  stand  was  then 
grown  to  age  100. 

Excellent  discussions  of  cover-forage  ratios,  hiding  cover,  and  thermal  cover 
requirements  for  deer  and  elk  are  presented  by  Thomas  and  others  (1979a). 
Thermal  cover  is  defined  to  be  any  stand  of  coniferous  trees  40  ft  or  more  tall, 
with  an  average  canopy  closure  exceeding  70  percent.  Figure  17  compares 
canopy  closure  and  tree  height  for  the  two  prescriptions.  The  shelterwood  stand 
reaches  a  top  height  of  56  ft  and  79  percent  canopy  closure  at  age  60.  Beyond 
60  years,  the  canopy  is  nearly  fully  closed,  providing  thermal  cover  throughout 
the  rest  of  the  projection.  The  clearcut  stand  supplies  less  adequate  thermal 
cover,  at  an  older  age  (80  years  and  beyond). 
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PERCENTAGE  CANOPY  CLOSURE  -  SHEUTERWOOD  AND  CUARCUT 
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MEAN  TREE  HEIGHT  (FEET)  -  SHELTERWOOD  AND  CLEARCUT  TREATMENTS 
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Figure  17.  — Thermal  cover  compared  tor  shelterwood  and  clearcut 
prescriptions  for  the  wildlife  example  stand:  (a)  stand  canopy  closure 
versus  stand  age  following  initial  thinning  in  1984:  (b)  average  tree 
height  versus  stand  age. 
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Hiding  cover  (defined  by  Thomas  and  others  [1979b]  to  be  vegetation  capable 
of  hiding  90  percent  of  a  standing  adult  elk  from  view  at  a  distance  equal  to  or 
less  than  200  ft)  can  be  compared  between  alternative  treatments  by  looking  at 
shrub  cover,  crown  profile  area  in  the  first  10  ft  of  height  in  the  stand,  the  sum 
of  stem  diameters,  and  trees  per  acre  (fig.  18).  Predicted  shrub  cover  is  44  per- 
cent for  the  shelterwood  treatment  and  43  percent  for  the  clearcut  treatment 
following  the  1984  thinning.  It  increases  to  52  percent  at  age  30  in  the  clearcut 
stand.  In  the  shelterwood  stand,  shrub  density  decreases  to  31  percent  at  age 
20,  and  then  the  overstory  removal  triggers  a  second  wave  of  increasing  shrub 
cover.  By  age  40,  shrub  density  is  again  equal  in  the  two  treatments,  at  about 
45  percent  cover. 

In  general,  hiding  cover  in  stems  and  tree  crowns  is  greater  in  the  shelter- 
wood  treatment,  primarily  because  regeneration  is  more  successful.  In  both 
treatments,  area  in  crowns  in  the  lower  10  ft  of  the  stand  begins  at  nothing  at 
stand  age  0  and  increases  to  a  maximum  at  about  30  or  40  years.  After  40 
years,  the  lower  canopy  level  begins  to  grow  above  the  height  where  it  can  be 
considered  effective  hiding  cover.  Beyond  60  years,  stem  area  contributes  more 
to  hiding  cover,  and  crown  profile  area  contributes  less. 

Figure  19  illustrates  how  canopy  development  predicted  by  COVER  may  be 
displayed  graphically  through  time.  The  vertical  distribution  of  crown  profile 
area  is  shown  for  the  two  prescriptions  immediately  before  thinning  in  1984, 
and  at  several  points  in  time  following  harvest.  Similar  stand  profiles  could  be 
drawn  to  represent  numbers  of  trees,  percentage  of  canopy  closure,  crown  vol- 
ume, or  foliage  density  by  height.  These  values  may  be  useful  in  relating  bird 
habitats  to  the  structure  of  vegetation  (for  example,  see  Langelier  and  Garton 
in  press  b;  Peterson  1982)  or  the  "life  form"  association  with  stand  successional 
stage  proposed  by  Thomas  and  others  (1979b). 

The  wildlife  example  presented  here  shows  how  stem  area,  crown  profile  area, 
and  related  values  can  be  interpreted  as  thermal  and  security  cover  for  big 
game.  Foliage-height  profiles  and  crown  volumes  by  height  classes  are  also  use- 
ful statistics  in  analyzing  bird  habitat  relationships.  The  interpretation  of 
canopy  and  shrub  height  and  density  into  wildlife  cover  is  hypothetical  for  the 
example  presented  and,  of  course,  depends  on  knowledge  of  actual  stand 
conditions. 

Vegetation  management  for  wildlife  use  requires  interdisciplinary  knowledge 
of  the  interactions  of  vegetation,  site  and  topographic  conditions,  silvicultural 
options,  and  road  and  harvest  operations.  The  combined  Prognosis/COVER 
model  can  provide  information  about  the  condition  of  the  vegetation,  including 
species  composition,  size,  and  distribution  of  both  the  overstory  and  under- 
story.  Thomas  (1979)  points  out  that  habitat  use  does  not  follow  some  arbitrary 
step  function,  but  that  wildlife  species  use  vegetation  despite  what  wildlife  biol- 
ogists define  to  be  less  than  optimum  conditions  (70  percent  canopy  closure  for 
thermal  cover,  for  example).  COVER  values  are  expressed  as  continuous 
through  time,  rather  than  as  threshold  values,  so  that  the  user  may  evaluate 
their  implications  for  wildlife  habitat  management.  COVER  output  can  be  writ- 
ten to  a  disk  file  (by  using  the  second  parameter  on  the  COVER  keyword)  for 
later  summarization  and  graphical  display. 
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Figure  18. —Hiding  cover  compared  for  shelterwood  and  clearcut  prescriptions  for  ttie  wildlife  example 
stand:  (a)  percentage  shrub  cover  versus  stand  age  following  initial  thinning  in  1984:  (b)  crown  profile  area  in 
the  lower  10  ft  of  the  stand  versus  stand  age:  (c)  sum  of  stem  diameters  versus  stand  age:  (d)  number  of 
trees  versus  stand  age. 


42 


SUM  OF  STEM  DIAMETERS  (HT)  -  SHELTERWOOD  AND  CLEARCUT  TREATMENTS 
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Figure  18.— (Con.) 
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Figure  19.  — The  development  through  time  of  crown  profile 
area  by  stand  height  for  the  wildlife  habitat  example. 
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Hydrologic 
Applications 


Forest  Insect  Pest 

Modeling 

Applications 


Succession  Modeling 
and  Planned 
Improvements 


The  COVER  extension  can  be  used  to  link  natural  or  human-made  vegetation 
changes  to  watershed  impacts,  primarily  through  predicted  canopy  closure  and 
shrub  cover  values.  The  sum  of  cover  conditions  across  stands  in  a  watershed 
can  be  expressed  as  the  percentage  of  ground  exposed  to  precipitation  and  run- 
off following  management. 

The  COVER  relationships  for  conifer  foliage  biomass  and  crown  shape  are 
currently  being  used  in  the  Western  Spruce  Budworm  Outbreak  Model 
(Sheehan  and  others  in  preparation).  The  foliage  equations  predict  total  poten- 
tial foliage  biomass  on  undamaged  trees.  Assumptions  about  the  partitioning  of 
foliage  by  age  classes  within  branches,  the  distribution  of  foliage  within  crowns, 
and  the  predicted  effects  of  defoliation  on  future  foliage  production  are  all  com- 
ponents of  the  western  spruce  budworm  (WSBW)  model.  These  relationships 
help  predict  insect  damage  to  trees  through  availability  of  food  and  pattern  of 
larvae  dispersal  vertically  in  the  stand. 

The  combination  of  understory  development  and  vertical  and  horizontal 
canopy  development  relationships  comprising  the  COVER  extension  represents 
a  framework  for  which  quite  detailed  successional  trends  can  be  displayed 
through  time.  Planned  future  studies  will  more  explicitly  link  the  shrubs  com- 
ponent to  the  regeneration  system  (Ferguson  and  others  1985)  and  small  tree 
development  models  (Wykoff  1985)  by  modeling  effects  of  shrub  competition  on 
small  conifer  establishment  and  growth  rates. 

A  second  improvement  planned  for  COVER  includes  the  option  of  making 
predictions  on  individual  sample  points  within  a  stand,  thus  allowing  a  heter- 
ogeneous site  to  be  represented  in  greater  resolution.  This  will  improve  the 
prediction  of  shrub  conditions,  allow  reporting  of  within-stand  variance  statis- 
tics, and  provide  a  measure  of  the  spatial  distribution  of  overstory  and  under- 
story cover. 

Work  is  progressing  on  a  graphical  display  link  that  shows  the  vertical  and 
horizontal  relationships  of  trees  and  shrubs  through  time.  The  display  is  in  the 
form  of  a  "lollipop"  diagram  in  which  each  of  several  figures  of  a  certain  height 
and  shape  represent  different  types  of  tree  and  shrub  records  in  the  stand. 

The  Prognosis/COVER  program  incorporates  models  that  are  specific  to  cer- 
tain species  and  conditions  prevalent  in  the  Northern  Rocky  Mountains,  but  it 
is  also  a  general  system  that  can  be  calibrated  to  local  conditions.  As  shrub 
and  crown  data  specific  to  other  areas  and  habitat  types  become  available,  new 
relationships  can  be  incorporated  into  the  model  to  expand  the  range  of 
predictions. 
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APPENDIX:  SUMMARY  OF  COVER  KEYWORDS 


Category  Keyword 


Keyword  use  and 
associated  parameters 


Default  parameter 
or  conditions 


Program 

control 


CANOPY 


COVER 


END 

SHRUBS 


Calibration  SHRBLAYR 


SHRUBHT 


SHRUBPC 


Output  DATELIST 
control 

DEBUG 


NOCOVOUT 

NOSHBOUT 
NOSUMOUT 

SHOWSHRB 


Compute  canopy  cover  statistics. 

Invoke  the  COVER  extension;  always 
the  first  keyword  record. 

field  1:  Cycle  in  which  COVER 

calculations  begin. 

field  2:  Dataset  reference 

number  for  COVER  output. 
Last  keyword  record;  return  control 
to  main  program. 
Compute  shrub  statistics. 

field  1:  Number  of  years  since 

stand  disturbance. 

field  2:  Habitat  type  code. 

field  3:  Physiographic  type  code, 
field  4:  Disturbance  type  code. 

Enter  shrub  calibration  values  by 
shrub  layer. 

fields  1,3,5:  Average  height  (ft) 

of  three  shrub  layers. 

fields  2,4,6:  Average  cover  (%) 

of  three  shrub  layers. 
Enter  shrub  height  calibration 
values  by  species. 

Up  to  four  supplemental  records: 

Species  code  and  height  (ft)  in 

consecutive  10-column  fields. 

Enter  shrub  cover  calibration  values 

by  species. 

Up  to  four  supplemental  records: 
Species  code  and  cover  (%)  in 
consecutive  10-column  fields. 

Print  date  of  last  revision  for  COVER 

model  subprograms  and  common 

areas. 

Request  printout  of  COVER 
calculations  for  tree  and 
shrub  list. 

field  1:  Cycle  in  which  debug 

output  is  to  be  printed. 
Suppress  the  canopy  cover  statistics 
display. 

Suppress  the  shrub  statistics  display. 
Suppress  the  canopy  and  shrubs 
summary  display. 

Select  additional  shrub  species  to  be 
displayed. 

Supplemental  record:  Species 
codes  right  justified  in  six 
consecutive  10-column  fields. 


None 
1 

18 


Stand  age  (STDINFO  card), 
or  3  years 

Stand  habitat  type  code 
(STDINFO  card) 
2  (Lower  slope) 
1  (None) 

No  calibration 


No  calibration 


No  calibration 


None 


Print  in  all  cycles 

Display  printed 

Display  printed 

Display  printed 
Print  nine  species 
which  account  for 
most  cover. 
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The  COVER  model  predicts  vertical  and  horizontal  tree  canopy  closure,  tree 
foliage  biomass,  and  the  probability  of  occurrence,  height,  and  cover  of  shrubs 
in  forest  stands.  This  paper  documents  use  of  the  COVER  program,  an  adjunct 
to  the  Stand  Prognosis  Model.  Preparation  of  input,  interpretation  of  output,  pro- 
gram control,  model  characteristics,  and  example  applications  are  described. 


KEYWORDS:  stand  structure,  crown  width,  crown  shape,  canopy  closure,  foliage 
biomass,  shrub  cover,  shrub  height,  shrub  occurrence,  stand 
simulation 


The  Intermountain  Research  Station,  headquartered  in  Ogden, 
Utah,  is  one  of  eight  Forest  Service  Research  stations  charged  with 
providing  scientific  knowledge  to  help  resource  managers  meet  human 
needs  and  protect  forest  and  range  ecosystems. 

The  Intermountain  Station's  primary  area  includes  Montana,  Idaho, 
Utah,  Nevada,  and  western  Wyoming.  About  231  million  acres,  or  85 
percent,  of  the  land  area  in  the  Station  territory  are  classified  as 
forest  and  rangeland.  These  lands  include  grasslands,  deserts, 
shrublands,  alpine  areas,  and  well-stocked  forests.  They  supply  fiber 
for  forest  industries;  minerals  for  energy  and  industrial  development; 
and  water  for  domestic  and  industrial  consumption.  They  also  provide 
recreation  opportunities  for  millions  of  visitors  each  year. 

Several  Station  research  units  work  in  additional  western  States,  or 
have  missions  that  are  national  in  scope. 

Field  programs  and  research  work  units  of  the  Station  are  maintain- 
ed in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana  State 
University) 

Logan,  Utah  (in  cooperation  with  Utah  State  University) 

Missoula,  Montana  (in  cooperation  with  the  University  of 
Montana) 

Moscow,  Idaho  (in  cooperation  with  the  University  of  Idaho) 
Ogden,  Utah 

Provo,  Utah  (in  cooperation  with  Brigham  Young  University) 
Reno,  Nevada  (in  cooperation  with  the  University  of  Nevada) 


